Quantum entanglement of a large number of photons 
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A bipartite multiphoton entangled state is created through stimulated parametric down-conversion 
of strong laser pulses in a nonlinear crystal. It is shown how detectors that do not resolve photon 
number can be used to analyze such multiphoton states. Entanglement of up to 12 photons is 
detected using both the positivity of the partially transposed density matrix and a newly derived 
criteria. Furthermore, evidence is provided for entanglement of up to 100 photons. The multi- 
particle quantum state is such that even in the case of an overall photon collection and detection 
efficiency as low as a few percent, entanglement remains and can be detected. 

PACS numbers: 42.65.Lm, 42.50.Dv, 03.65.Ud 



In recent years, small numbers of entangled particles 
have been used for testing quantum mechanics and for 
implementing various quantum information protocols I']. 
However, other tests probing the validity of quantum 
decoherence models and additional quantum infor- 
mation protocols will require entangled states of large 
numbers of particles. Bipartite multiphoton states, the 
subject of this paper, can be used to test foundations of 
quantum theory[3,0|j and for quantum cryptography'^. 
Furthermore, it has been shown that phase sensitive 
measurements!^ and quantum photolithography0 can 
exceed classical boundaries imposed by the wavelength 
of hght, by using multiple entangled photonsQ. 

In this Letter we demonstrate the generation of a 
bipartite entangled state of many photons in two spa- 
tial modes, as produced by stimulated parametric down- 
conversion (PDC)|i,[ll[llII2- The Hamiltonian for the 
generation of polarization entangled photons |^ is given 

by 
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Horizontally (h) and vertically (v) polarized photons oc- 
cupy two spatial modes (a and &); k is a coupling con- 
stant that depends on the nonlinearity of the crystal and 
the intensity of the pump pulse. The resulting photon 
state is given bvjllj 
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where, for example, \m)a^ represents m vertically polar- 
ized photons in mode a. The interaction parameter r 
depends linearly on the crystal length and on k. The 
state IV') is a superposition of the states |V'n ) of i^" 
distinguishable photon pairs. Each is an analog of 
a singlet state of two spin-n/2 particles, thus l^p) is in- 
variant under joint rotations of the polarization bases 



of both modes. The average photon pair number is 
= 2sinh^T. Although photons are created in pairs, 
the resulting state cannot be factorized into individual 
pairs. As a result of the stimulated emission process, the 
pairs are indistinguishable such that they form a single 
multiphoton entangled state. Previous PDC experiments 
have been restricted to t < 0.1, resulting in the detection 
of at most 4 to 5 photons in an entangled state [l3.ll^ll3 |. 
This work addresses the region of r > 1, where entangled 
states of large numbers of photons can be generated. 

Our setup is switchable between single-pass and 
double-pass jl3 of a pump pulse through a BBO nonlin- 
ear crystal (see Fig. ^ . The pump is a frequency-doubled 
amplified Tiisapphire laser, giving 200 fs pulses with 5 /xJ 
per pulse. Two down-converted modes a and b are se- 
lected by coupling into single-mode fibers through 5 nm 
narrow bandpass filters, justifying the use of the Hamil- 
tonian of Eq. n For each spatial mode, two orthog- 
onal polarizations are separated by a fiber polarization 
beam splitter and detected by silicon avalanche photodi- 
odes (APD). Single-photon count rates and coincidence 
rates were recorded as functions of the pump power in 
three polarization bases: horizontal/vertical linear (hv), 
plus/minus 45° linear (pm) and right/left circular (rl). 

To characterize states of multiple photons, it is desir- 
able to have detectors that can resolve photon number. 
Although such detectors exist [Hill, Hill Till, the pho- 
ton number resolution is always limited by losses. There- 
fore, to determine the actual multiphoton state produced 
at the source, it is necessary to perform a probability 
analysis of the experimental data, based on the physics 
of the detection scheme. 

We use APDs that give no direct information about 
the number of detected photons. Nevertheless, the 
probability to obtain a signal depends on the photon 
number 20, 21j. For a total collection efficiency rj (a com- 
bination of the APD detectors and the optical coupling 
efficiencies) , the triggering probability given an m-photon 
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FIG. 1; The experimental setup. The pump pulses pass twice 
through the nonlinear BBO crystal, with a controlled delay 
between the passes. Down-converted photons from the first 
pass are re-injected into the crystal together with the pump 
second pass 12]. The configuration can be switched to single- 
pass by blocking the path of the first down-conversion. Addi- 
tional BBO crystals are inserted in order to compensate for 
temporal walk-off. The two spatial modes are collected into 
single-mode fibers through narrow bandpass filters. Photons 
in the different modes are detected by four APDs. 



state is 



1 - Po = 1 - (1 - 77)^ 



(3) 



where Pq is the probability of not detecting a photon. 

For the state of Eq. [3 the detection probabihties for 
a single spatial mode and coincidences between any two 
modes as functions of r and rj were derived. For example, 
the probability per pulse to trigger a single detector is 
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We measured probabilities as functions of the pump in- 
tensity / up to a maximal intensity Imax, and fitted the 
results with the collection efficiencies rji of the four modes 
ah, ciy, bh and by, and the maximum interaction param- 
eter Tmax defined as Tmax = T^lmaxll- 

Single-photon counts of one polarization mode and co- 
incidence counts are presented in Fig. This data is 
from a single-pass configuration. All the results in various 
polarization bases were successfully fitted with the same 
parameters Tmax = 2.30 ±0.05 and = 1.9 ±0.2% for all 
four modes, strongly supporting the model of Eq. Q. 
The stimulated emission enabled the direct observation 
of coincidences that can only occur from events of at least 
three or four pairs. By collecting only one polarization 
and splitting the photons into two detectors with a beam 
splitter, we counted coincidence events of the form ah-au- 
bh and of the form ah-au-hh-bh that originated from three 
(or more) and four (or more) photon-pairs, respectively. 
Figure|2b combines all the measurements as a function of 
T. The larger the number of relevant pairs for an event. 
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FIG. 2: (a) Experimentally observed single photon count 
rates of the mode (squares), one-pair coincidence rates 
between ah and bv (circles) and two-pair coincidence rates 
between and bh (triangles) as a function of the pump 
pulse energy. The one-pair coincidences can arise from one 
or more pairs, while the two-pair coincidences can arise from 
two or more. Fits are included as solid lines, (b) The r de- 
pendance of the results of Fig. and coincidence events that 
can be generated by at least three pairs (inverted triangles) 
and four pairs (diamonds) of photons. 



the steeper is its graph, as expected from a multiphoton 
stimulated process. The slopes for t < 1 of the single 
photon counts and the one-pair coincidence are parallel 
and linear with pump pulse energy, as both result from 
one-pair events. All the graphs should saturate for large 
enough r at the repetition rate of 20 kHz. The maxi- 
mum interaction parameter achieved corresponds to 100 
photons per pulse on average. 

We have shown a stimulated emission process in which 
many photons are created in a way consistent with the 
state of Eq. |21 This does not yet prove the specific quan- 
tum correlations between the photons described by that 
state. We will now present two criteria that prove the 
presence of entanglement. Our approach is to use a low 
overall detection efficiency such that in the relevant pa- 
rameter regime we detect at most two photons. We show 
that entanglement is still present in this situation. This 
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FIG. 3; Measured density matrices in the (1,1) subspace for 
r = 0.2 (a) and r = 1.85 (b). Only the real part of p is shown 
as its imaginary part is negligible. 
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proves entanglement not only of the detected photons, 
but also of the initial state before most of the photons 
were lost, as it is impossible to form an entangled state 
by applying a local operation to unentanglcd photons. 

The Positivity of the Partially Transposed (PPT) den- 
sity matrix is a separability (non-entanglement) crite- 
rion for bipartite systems, such as the PDC system stud- 
ied here|22|. Consider the density matrix p in the sub- 
space where only one photon after losses is detected in 
each of the spatial modes. This restriction will be justi- 
fied below. The total probability for such a detection is 
= Phv + Pvh + Phh + Pvv, where Phv is the proba- 
bility to detect a coincidence between an and 6^, etc. By 
only considering events from this (1,1) subspace, prob- 
abilities can be normalized as pij = Pij/P(i.i). We use 
\hh), \vh) and \vv) as the basis states for p. We 

also define the single-pair visibility V for different polar- 
ization detection bases. For example, in the hv basis for 
mode a and the pm basis for mode b 
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The elements of p can be readily obtained from 
combinations of visibilities, a process known as state 
tomography [2^. Density matrices, as measured for low 
(r = 0.2) and high (r — 1.85) values of t, are pre- 
sented in Fig. 13 The measured density matrices are 
consistent with the state of Eq. |2J for low tau the den- 
sity matrix approaches the familiar two-photon IV'r) ^^^^ 

state {phv,hv = Pvh,vh = —Phv,vh = —Pvh,hv — 1/2), for 

high tau hh and vv coincidences are detected {phh,hh and 
Pvv,vv are no longer small) as a result of multiple photon 
pairs before detection losses. Considering these domi- 
nant 6 terms, the partially-transposed matrix p^-^ will 
only have positive eigenvalues if 



16 • PhhPvv 



>1. (6) 



The violation of the above separability criterion proves 
entanglement. For example, for a pure {ipi) state Ci is 



FIG. 4: The experimentally measured entanglement criteria. 
To detect entanglement the PPT criterion Ci (circles) must 
be smaller than one, while the visibility/spin-correlation cri- 
terion C2 (squares) must be larger than one. Both criteria 
detect entanglement up to an interaction of r = 1.3 (dashed 
line), corresponding to a state with 12 photons on average. 



zero. The circles in Fig. 0]show the measured Ci as a 
function of r. Entanglement is detected up to an inter- 
action of r = 1.3, corresponding to 6 indistinguishable 
photon-pairs on average. For this r value and detection 
efficiency of 2%, the ratio of the probabilities for a detec- 
tion from a higher photon number subspace and from a 
(1,1) subspace is about 0.06. This low ratio justifies con- 
sidering only (1,1) subspace events for the density matrix. 

We now derive a second separability criterion tailored 
to detect the type of entanglement created by PDC. The 
visibility (Eq. [SJ can be rewritten as the spin anti- 
correlation between the two spatial modeSjlOj: 



Vhv,hv = Phv + Pvh - Phh - Pv 



(7) 



where cr^ are the Pauli matrices. The total spin correla- 
tion is: 
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A product state, is maximally correlated (anti- 
correlated) when the two spins are parallel (anti-parallel). 
It is convenient to rotate the two correlated spins to one 
of the principal bases. The correlation in that basis will 
be ±1, and zero in the other two bases. A general sepa- 
rable state (a mixture of product states) can not have a 
higher correlation. On the other hand, an entangled state 
such as the Bell \tpi) state can have all the correlations 
as -1. Thus, an upper bound criterion for a separable 
state in the (1,1) subspace is 
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FIG. 5: One-pair visibility as a function of the interaction 
parameter for the polarization bases hv (solid circles), pm 
(open circles), and rl (crossed circles), and for the model 
fit (solid line). The prediction for the ansatz state of distin- 
guishable entangled pairs (dashed line) represents the vis- 
ibility upper bound for uncorrelated pairs of photons. The 
observed results are above this bound, indicating the indis- 
tinguishable nature of the generated photon-pairs. 

The measured values of C2 are shown as squares in Fig. 
0] This second criterion detects entanglement also up to 
r = 1.3, or up to about 6 photon pairs. Although based 
on different arguments, the two criteria detect entangle- 
ment up to the same interaction values. This boundary is 
mainly set by the limitations of the APDs and not by the 
actual entanglement present in the generated state. One 
can show theoretically that some entanglement remains 
for arbitrarily high r and losses 24] . 

To provide experimental support for entanglement over 
the entire detected interaction range, the one-pair vis- 
ibilities in three polarization bases as measured with a 
double-pass setup, optimized for collection efficiencies of 
77 = 9 ± 0.7%, are shown in Fig. Oand compared to their 
theoretical prediction (solid line). The r-dependence 
is approximately the same in the different polarization 
bases, consistent with the state rotation invariance. For 
comparison, the visibility for an ansatz state of distin- 
guishable pairs of entangled photons was also calculated 
(see dashed line in Fig. [SJ. We used the same pair- 
number distribution as in PDC, but assumed that the 
different pairs occupy different modes. The predicted 
visibility for this ansatz state is considerably lower than 
the PDC visibility curve and the experimental results. 

In conclusion we have demonstrated the generation of a 
bipartite state of up to 50 indistinguishable photon pairs. 



Entanglement up to 12 photons has been proven while 
evidence has been given for entanglement up to 100 pho- 
tons. We have shown that it is possible to explore quan- 
tum entanglement even after the state suffered significant 
losses and even with detectors that have limited photon 
number resolution. The studied multiphoton entangled 
state is of particular interest for quantum cryptography, 
quantum metrology and for tests of the foundations of 
quantum mechanics with large-spin systems. 
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